Charge density plateaux and insulating phases in the t — J model with ladder geometry 
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We discuss the occurrence and the stability of charge density plateaux in ladder-like t — J systems 
(at zero magnetization M = 0) for the cases of 2- and 3-leg ladders. Starting from isolated rungs 
at zero leg coupling, we study the behaviour of plateaux-related phase transitions by means of first 
order perturbation theory and compare our results with Lanczos diagonalizations for t — J ladders 
(N = 2x8) with increasing leg couplings. Furthermore we discuss the regimes of rung and leg 
couplings that should be favoured for the appearance of the charge density plateaux. 
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I. INTRODUCTION 



The t — J model has been introduced as the first order 
correction of the extreme atomic limit of the Hubbard 
mode! 1,2 and is considered the simplest model including 
the low energy physics of doped ladder systems The 
phase diagram of the two leg Hubbard model has been 
investigated in |5j by means of a renormalization group 
approach valid for small values of the on-site Coulomb in- 
teraction U but for arbitrary charge density and arbitrary 
hopping along the rungs. The phase diagram is classified 
as different C x S y phases which denotes x gapless charge 
modes and y gapless spin modes. They have shown under 
which condition a phase of C±So appears which is analog 
of either a superconductor or charge-density wave. The 
extension of this approach to the iV-leg Hubbard model 
can be found in |6j where the dimensional crossover as 
N — ► oo is discussed. The charge and spin gap for the 
two leg Hubbard model and its dependence on the on- 
site Coulomb interaction and the rung hopping parame- 
ter t have been calculated by means of the density matrix 
renormalization group (DMRG) and compared with the 
previous weak-coupling RG. 7 The effect of an additional 
nearest neighbour Coulomb repulsion V has been stud- 
ied in A charge order (metal- insulator) transition 
was found at charge density p = 1/2 from a homoge- 
neous state to a charge density wave. The influence of 
an anisotropy between leg and rung couplings in Hub- 
bard and t — J models on specific correlations, which sig- 
nal the metal-insulator transition, has been investigated 
in 0. The metal-insulator transition is accompanied by 
the opening of a gap, which appears as a plateau in the 
charge density p{p) as a function of the chemical poten- 
tial \x. The charge density plateau in p{p) looks similar 
to the plateaux in the magnetization curve M(B) found 
in the spin ladder systems and one might ask whether 
the mathematical foundations for both plateaux are the 
same. This is indeed the case and becomes evident if 
one maps the t — J Hamiltonian on a spin-1 Hamilto- 
nian with broken SU(3) symmetry— The Lieb-Schultz- 
Mattis theoremii has been extended to quasi onedimen- 
sional fermionic systems 12 and the momenta of low-lying 
excitations could be classified thereby. The quantization 



rule of Oshikawa, Yamanaka and AfTled*^ predicts as 
well possible plateaux in the charge density p(ju). 

The values of the charge density (and magnetization) 
at the plateaux are fixed by the geometry of the system 
(e.g. the number of legs in a ladder system). 

In this paper we will study the ground states of the 
t — J Hamiltonian on a ladder with n legs (n = 2, 3) 



H M = fH ln] + t , H [n\ 



(1.1) 



t, t' are the hopping parameters, and 



N r 



#N = ^2hW(x,a), a = J/t, 



(1.2) 



N, — -1 



H l n] = J2 M nl (*»«') a 1 = j'/t' 



define the contributions of the couplings along the rungs 
and legs, respectively. The spin exchange is included in 
the ratios a = J/t, a' = J'/t'. The latter are depicted 
in Figs. and [21 for the cases of a two and a three leg 
ladder. 



N+2 N+3 



2N-2 2N-1 




FIG. 1: Structure of a two leg ladder with open boundary 
conditions and 2iV,. sites. 

A charge density plateau in p(p) - where p, is the chem- 
ical potential - is signalled by discontinuous changes in 
the slope of the ground state energy per site as function 
of p [cf . eq. H2.12fl ] . They emerge immediately in the "lo- 
cal rung approximation" 9 (or "bond operator theory"—) , 
where the ground states are direct products of rung clus- 
ter states. The corresponding energies turn out to be 
piecewise linear in p. 
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It is the purpose of this paper to go beyond the local 
rung approximation by means of a systematic perturba- 
tion theory in the leg coupling t' . To first order, this 
leads to an effective interaction between the rung cluster 
states. The corresponding effective Hamiltonians are de- 
fined on a chain with N r sites - where N r is the number 
of rungs - and can be diagonalized numerically. 

The outline of the paper is as follows: In Section Ull we 
first treat the two leg ladder. Here, the effective Hamil- 
tonians for p < 1/2 and p > 1/2 look like modified t — J 
models on a chain with N r sites. 

In Section ITTT1 we extend all considerations to the three 
leg ladder, where plateaux appear at p = 1/3 and p = 
2/3. First order perturbation theory in if leads to three 
different effective Hamiltonians in the regimes < p < 
1/3, 1/3 < p < 2/3, 2/3 < p < 1. 

In Section lTvl we compare the perturbative results with 
Lanczos diagonalizations on ladder systems and discuss 
consequences for the phase diagram. 



A. Oth order perturbation theory in t' 

The lowest energy eigenstates of the rungs (x, x + N r ) 
with Hamiltonian h(x, x+N r , a) can be easily calculated. 

\Q(x) = 0) = |0,0) 

|Q(aO = l) = -^(|a,0) + |0,a)) (2.6) 

\Q(x) = 2) = _)_|-,+» 

The right-hand side defines the states on the sites x 
and x + N r of the rung; means hole with charge zero. 
Q(x) = 1 means one electron with spin a = ±1. Q(x) = 2 
represents an electron pair coupled antiferromagnetically 
to a total spin 0. Therefore, we have charges 0, 1 and 2 
for \Q(x) — 0), \Q(x) = 1)) and \Q(x) = 2) respectively. 
The eigenvalues are given by the following equations. 

h(x,x + N r ,a)\Q{x)) = e n \Q{x)) (2.7) 



II. T-J MODEL ON A TWO LEG LADDER 

According to the notation given in Fig. ^ the cou- 
plings along the rungs and legs which enter into the t — J- 
Hamiltonian on a two leg ladder are given by: 



e — 0;ei — — V.E2 = -a 



(2.8) 



(x, a) 
hf\x,a') 



h(x,x + N r ,a), (2.1) 
h{x,x + l,a') (2.2) 
+h(x + N r , x + 1 + N r , a') . 



For our purposes it is convenient to represent these 
couplings in terms of "constrained" permutation opera- 
tors 



h(x,y,a) = -P i(x,y) + 



in 



(x,y) (2.3) 



Here, Po\{x,y) is a permutation operator, which per- 
mutes the states at x and y, if they are occupied by 
one hole (0) and one electron (1). This operator allows 
the hopping of electrons and holes and forbids the double 
occupancy of each site x with two electrons. 

If both sites x and y are occupied with electrons, the 
operator 

Qn(x,y) - P(x,y)-l, (2.4) 
Qn(x,y)\x,y) = \y,x)-\x,y) (2.5) 

first interchanges the electrons at x and y; afterwards the 
original state is subtracted. 



In the limit of vanishing leg couplings t' = 0, the system 
of decoupled rungs has the following ground state. 



(1) For p = JjL - 



and 



I 2_ < I- o = 012 

2N r 2 N r — 2 ' y > ' > 1 

a < 2 the ground state is a direct product of (N r — 2q) 
rung states \Q(x) — l(cr)) and 2q states with \Q(x) = 0). 

Note that for a < 2, the creation of a pair state 
|Q(:r) = 2) and a charge zero state \Q(x) = 0) from 
two charge 1 states \Q(x) = l(cr = +1)) and \Q(x) — 
l(er = —1)) is not energetically favorable, since 

e +£2>2ei (2.9) 

Therefore, the ground state energy turns out to be 

E^ ] {p,a) = {N r -2q)e 1 + 2qe 

= -2N r p (2.10) 

(2) For p = 2§- = i + £ > |; q = 0, 1, 2, . . . and 
a < 2 the ground state is a direct product of (N r — 2q) 
rung state \Q(x) = 1) and 2q states \Q(x) = 2) with 
electron pairs. In this case the ground state energy is 

E^ ] (p,a) = (Nr-2q)e 1 +2qe 2 

= 2N r p(l - a) - (2 - a)N r 

(2.11) 

Note, in both regimes p < 1/2 and p > 1/2 the ground 
state energy is linear in p. The first derivative of the 
ground state energy 



d E 
d^2AV 



(2.12) 
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is related to the chemical potential which has a disconti- 
nuity at p = i . 

1 — a for p > 5 , v ' 

This is the first indication of a charge density plateau at 
p = 1/2 in the t — J model on a two leg ladder. 

To zeroth order in the leg coupling (i'), the eigenstates 
of the t- J Hamiltonian are product states of N r rung 
states: 

|Q(1), Q(2), . . . , Q(N r )) = J] |Q(as)) (2.14) 

X 

The rung quantum numbers Q(x) = 0, 1(<t = ±1), 2 are 
subjected to the conservation of total charge Q, 

Y,Q{*)=Q (2-15) 

X 

Moreover, the total spin, which originates from the spin 
of charge 1 states \Q{x) = 1), is assumed to be zero here: 

J2^q(x),i=0 (2.16) 

X 

Therefore the ground states with energies l|2.10|l . (|2.11|) 
are highly degenerate for p ^ 0, 1/2, 1. 

B. 1st order perturbation theory in t' 

A first order perturbation theory demands a computa- 
tion of the transition matrix elements: 

(Q'(l), Q'(2), . . . , Q'{N)\Hf ] |Q(1), Q(2), . . . , Q{N)) = 

N r -1 

A(Q'{x), Q'(x + l)\hf\x, a')\Q{x),Q{x + 1)) 

x=l 

(2.17) 

where 

A = II S Q'(v),Q(y)' 

yjtx,x+l 

and the diagonalization of the resulting effective Hamil- 
tonian on a chain with iV r -sites. The matrix elements of 
the leg operators between the rung states Eq. Ij2.6|l 

(Q'(l),Q'(2)\hf\l,a')\Q(l),Q(2)) = 

(Q'(1),Q'(2);Q(1),Q(2)) , (2.18) 

- the explicit spin dependence of the electrons has been 
omitted here - are listed in the following equations: 

(0,0; 0,0) = 0, ; (a, 0; a, 0) = , (2.19) 
(^^(Ti.o-a) = ^(ol,o 2 |P(l,2)-l|o 1 ,o 2 ),(2.20) 

(o,0;0,o) = -1, ; (a, 2; 2, a) = -~ , (2.21) 

(2, o- 2, a) = ~, ; (2, 2; 2, 2) = ~ (2.22) 



From these equations we can read off the effective Hamil- 
tonian on a chain with N r sites. Since the ground state 
structure differs for p < h and p > h as described in 
Sec. (II), we have to consider these two cases separately. 

(1) For p < 1/2 only rung states with Q(x) = and 
Q(x) = 1 are involved. From Eqs. (|2.19(l and l|2.20|l we 
see that the effective Hamiltonian 

N r -i 

ffi(a') = J2 h{x,x + l,a') (2.23) 

x=l 

with couplings: 

a' 

hi(x,x+l,a') = -P Q1 (x,x + l) + —Q 11 (x,x + l) (2.24) 

can be identified with a i — J model l|2.3|l on a chain with 
coupling a'/2, which is just half of the leg coupling. Note 
also that the charge density on the chain with N r sites 

Pi = §- = 2p (2.25) 

is just twice the charge density of the two leg ladder. 
In first order perturbation theory we get for the ground 
state energy on the two leg ladder: 

E (p,t — l,a;t \a) = —2N r p 

+t'E 1 { Pl = 2p, ai = (J j2) 
+0{t' 2 ) (2.26) 

Here £a(pi, a\) is the ground state energy of the "effec- 
tive" t — J Hamiltonian 12.23fl on a chain with N r sites. 
According to Eq. 1)2. 12|) we can calculate the chemical po- 
tential from the first derivative with respect to the charge 
density: 

pi(p, a; t\ a 1 ) = -1 + t' pi(pi, «i) (2.27) 

where 

. , . 1 dE\ , , 

^ Pl,ai) = W r ~dpl ( } 

is the chemical potential of the t — J model on a chain 
with N r sites and ot\ = a' /2 at p\ = 2p. 

(2) For p > 1/2 only rung states with Q(x) = 1(<t = 
±1) and Q(x) = 2 are involved. From Eqs. (|23U|) - (|2~2"2l) 
we see that the first part of the effective Hamiltonian: 

N r -1 

H 2 (a') = h 2 (x,x + l,a') (2.29) 

x=\ 

with couplings: 

a' 

h 2 (x,x + I, a) = {— P 2 i(x,x+ 1) + — Qn(x,x + 1) 
+D 2 (x,x + 1)} (2.30) 
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is indeed a t — J model, if we treat electron pair states 
\Q(x) — 2) as quasi-holes. The third term (D 2 ) on the 
right-hand side of l|2.30jl : 

(Q(x),Q(x + l)\h 2 (x,x + l)\Q(x),Q(x + l)) 



2N+1 2N.+2 2N+3 3N-2 3N-1 3N„ 






Q(x) 


a 
2 


Q(x) 




Q(x) 



(2.31) 

takes into account the non-vanishing diagonal terms 
(12.22(1 . The ground state energy E 2 (p 2 ,a 2 ) of the ef- 
fective Hamiltonian 



H 2 = ^ h 2 (x, x + 1, a') 



(2.32) 



on a chain of N r sites fixes the first order perturbation 
correction to the ground state energy of the two leg ladder 
system: 

E(p,t = l,a;t',a') = -2N r [l - a/2 - p(l - a)] 



t -E 2 (P2,a 2 ) + 0(t' 2 ) 



where 



P2 = 2(1 - p) , a 2 = a' . 



(2.33) 



(2.34) 



Finally, we get the following relation between the chem- 
ical potential p(p, a; t 1 , a') of the two leg ladder 



ft 2 (p 2 ,a 2 ) = 



1 dE 2 
N r dp 2 



(2.35) 



and for the effective Hamiltonian i|2.32[l on a chain of N r 
sites: 

t' 

p 2 {p, a;t,a) = 1 - a - —p 2 {p 2l a 2 ). (2.36) 

Combining l|2.27|l and l|2.36|l . we get for the width of 
the charge density plateau at p — ^ in the first order 
perturbation theory: 

W(a,t\a) = p 2 - pi 

= 2-a-t'\p, 2 {p 2 = l,a')/2 

+p x { Pl = l,a'/2)] (2.37) 



III. CHARGE DENSITY PLATEAUX ON A 
THREE LEG LADDER 

The perturbation treatment of the leg couplings will 
be applied now on the three leg ladders. The geometry 
and the notion of states can be seen in Fig.©. 



t,a 
N r +1 
t,a 



t ,ot 









2N. 



N-2 N-l N, 



FIG. 2: Structure of a three leg ladder with open boundary 
conditions and 3iV r sites. 



The Hamiltonian (|l.l|l for the 3-leg (n = 3) case is 
again constructed from t — J couplings l|2.3|l on the rungs 
and the legs, respectively, and is given by the components 

(x, a) — h(x, x + N r , a) + h(x + N r , x + 2iV,., a) , 
3 



hf\x,a') = ^h{x + {l-l)N r ,x + l + (l-l)N r ,a') 
1=1 

(3.1) 



A. Oth order perturbation theory in t' 

The lowest energy eigenstates on the rung Hamiltonian 

[31 1—1 

h r (x, a) have been calculated in |15j . 
\Q(x)=0) = |0,0,0) 

\Q(x) = l) = i{|a,0,0)+V2|0 ;( T,0) + |0,0,a)} 

\Q(p) = 2) = -j=L= {\a,+,-)-\o,-,+) 

+6|+,0,-)-6|- 0,+) 
+1+,-, 0)-|- +,0» 

\Q(x)=3) = -j={\a,a,-a) - 2\a,-a,a) 

+ | -a, a, <t)} 



with 



b = — = (Va 2 +8-a)/2. 

£ 2 

a = ±1. 



(3.2) 
(3.3) 



We have four types of states with charges Q(x) — 0, 1, 2, 3 
respectively. Even charge states with Q(x) = 0, 2 carry 
total spin and look like "composite bosons". Odd 
charge states with Q(x) = 1,3 have total spin 1/2 and 
look like "composite fermions" . The energies of the four 
states are: 

£0 = 0, £1 = -V2 

£ 2 = --(Va 2 + 8 + a), e 3 = --a. (3.4) 
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If we compute the ground state energy of the three leg 
ladder in the limit of vanishing leg coupling if = 0, we 
have to discriminate the following three cases: 

(1) regime: < p < 1/3 (p = 1/3 - 2q/3N r ) 

The ground state is a direct product of N r — 2q rung 

states \Q(x) — 1(ct)) with charge 1 and 2q rung states 

\Q(x) = 0) with charge 0. 

E[ 0) (p, a) = (N r - 2q)ei + 2qe Q = 3N r p £l (3.5) 

The creation of a charge 2 state (\Q(x) — 2) ) and a charge 
state (\Q(x) = 0)) from two charge 1 states costs en- 
ergy, since 



e + £2 > 2e\ for a < —= (3.6) 

v2 

(2) regime: 1/3 < p < 2/3 (p = 1/3 + 2q/3N r ) 
The ground state is a direct product of N r — 2q rung 
states with charge 1 and 2q rung states with charge 2. 

^ 0) (p,a) = (JV r - 2q)e l + 2qe 2 

= N r {(2- 3p)s 1 + (3p-l)e 2 } (3.7) 

The creation of a charge 3 state and a charge 1 state from 
two charge 2 states costs energy since 



£1 + £3 > 2e 2 for a > 



(3.8) 



(3) regime: 2/3 < p < 1 (p = 1 - 2q/3N r ) 
The ground state is a product of 7V r — 2q rung states 
with charge 3 and 2q with charge 2. 

£3°W) = (^r - 2g)e 3 + 2<7£ 2 (a) 

= iV r {(3p-2)e 3 + 3(l-p)e 2 } (3.9) 

Note that the ground state energies are linear again in 
the charge density. Therefore, we get for the chemical 
potentials. 



(3.10) 



p(p) 


1 dE 




3N r dp 




pf\p, a) 


= £1 


for 


P2 \P, a ) 


= £2-£\ 


for 


P^ (P= °) 


= £3 -£2 


for 



These are the results for the zeroth order (t 1 = 0) in the 
leg couplings [hi in i|3.1[) ]. 



B. 



1st order perturbation theory in t' 



The first order calculation starts from the matrix ele- 
ments of the leg couplings l|3.1|l between the zeroth order 





r 


< p < i 

— > — 3 


k <P< 1 

3 — r 1 — 3 


2 < p < 1 

3 _ p _ 


Cyni n O" 


Q = 0, 1 


Q = 1,2 


Q = 2, 3 


f^Oth 


£1 


£2 — £l 


£3 — £2 


Peff. 


pi = 3p 


P2 = 2 - 3/9 


pi = 3p - 2 


H eS. 


ii-ffi(ai) 




£3^3(03) 


eu. 


tx = t' 


4(^+2) 


+ . _ 3 +/ 

2(^+2) 


a eff. 


Qi = fa' 


3(6 2 +2) / 
2 2(b 2 +2^6+2) 


2(b 2 +2) / 

03 = 3 a 


plst 


tipi(pi,atx) 


— tl{l2{p2, OI2) 


tzP%(.P3,a.%) 



TABLE I: Effective couplings for a three leg ladder 



states (|3.2|l . The results of this tedious calculation are 
summarized in Table [I] 

The first row defines the 3 regimes of charge density p. 
In each regime, the zeroth order ground state is built up 
from direct products of rung cluster states with charge Q 
listed in the second row. The chemical potentials (|3.10(l - 
(|3.11|) in zeroth order are listed in the third row. 

First order perturbation theory in the leg couplings 
leads to the effective Hamiltonians with nearest neigh- 
bour interactions on a chain of N r sites in each sector of 
p (5-th row of Table P). These Hamiltonians contain two 
parts: 

N r -1 

tjHjiaj) = tj ^2 hj(x,x + l,aj), (j = 1,2,3), 

x=l 

= tj[H t _j(t = l,o = a 3 ) + Dj] . (3.12) 

The first one is of the t — J type l|2.3(l with effective hop- 
ping parameter tj (6-th row of Table [IJ and spin coupling 
aj (7-th row of Table HJ. 

The second part Dj,j = 1,2,3 with D\ = takes into 
account diagonal terms, which are not present in the t — J 
Hamiltonian: 



{l,2\h 2 {x,x + l,a')\l,2) 



-(^+3) 
8(6 2 + 2) 



(3.13) 



-(h 4 + 2h 2 + 3) 
(2,2\h 2 (x,x + l,a')\2,2) = ■ ^ V (3.14) 



2(6 2 + 2) 2 



{2,3\h 3 {x,x + l,a')\2,3) = 
(3,3\h 3 (x,x+l,a')\3,3) = 



-a' 
~2~ 
-_of_ 
~2~ 



(3.15) 
(3.16) 



Note that in the three leg ladder case the effective hop- 
ping terms i 2 and t% as well as the effective couplings a 2 
and «3 depend on the rung coupling a via l|3.3l) . This 
does not occur in the two leg ladder case. 

If we denote the ground state energies of the effective 
Hamiltonian Hj on a chain with N r sites by Ej(flj,ctj) 
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and the corresponding chemical potential by 

1 dEj 



(3.17) 



we can express the first order correction to the chemical 
potential of the three leg ladder in terms of l|3.17[) (last 
row of Table P). The relation between the charge density 
p on the ladder system and the charge density p 3 in the 
effective one dimensional system can be found in the 4-th 
row of Table D 



IV. NUMERICAL RESULTS 

In this section we are going to present numerical results 
for the ground state energies and the chemical potentials 
of the two and three leg ladder. Our results were ob- 
tained with open boundary conditions to facilitate the 
comparison with future DMRG calculations which can 
be done on larger systems. Other boundary conditions - 
e.g. periodic ones - can be incorporated as well. 



A. Two leg ladders 

The ground state energies Ej(pj,aj) of the effective 
Hamiltonians Hj j = 1,2 [cf. (j2~2"3)) and l(2~2"5)l ] on 

a chain with N r sites have been computed for N r — 
8, 10, 12, 14, 16, (18) 1 and 

a' 

a = 2.7, ai = — = 1.35, a 2 = a' = 2.7 . (4.1) 

For p\ = and P2 — these energies are known: 

E x {pi = Q, ai ) = 0, (4.2) 



E 2 (p 2 = 0,a 2 ) = -a 2 (N r - 1) 



(4.3) 



For pi = 1 and p 2 = 1 the ground state energies are 
given by the nearest neighbour Heisenbcrg chain (with 
open boundary conditions) 

E 1 {p 1 = l,a 1 =a'/2) = ^E 2 (p 2 = 1, a 2 = a')(4A) 

The finite-size dependence has been analyzed with an 
ansatz: 

for the ground state energies per site. 



1 The 18-site systems have been evaluated for Q = 0, 2, 4, 6, 16, 18. 
For the Q values inbetween, the dimension of the Hilbert spaces 
exceeded our computer capacities. 



From H4.3fl we get 

N 2 ( Pj = 0) = -1 



(4.6) 



A finite-size analysis of the Heisenberg chain [t — J chain 
at p — 1) yields 



N 3 { Pj = 1) 



-0.6. 



J = 1,2, (4.7) 



which means that finite-size effects change with pj. We 
assume here a linear interpolation: 



-l + 0.4 Pi , j = 1,2, (4. 



between the boundary values (|4.6|l and (|4.7|l . 

Indeed, this procedure has the effect, that the data 
points for N — 12, 14, 16, (18) follow unique curves 
ei(pi,ai) and e 2 (p 2l a 2 )/2 as is demonstrated in Fig-EJ 
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FIG. 3: Ground state energies per site bji.j(pj,aj) [Eqs. 
114.511 and 14.9H 1 for the effective Hamiltonians Hj (j = 1, 2) 
of the two leg ladder with Q2 = a' — 2.7 - results for 
N r = 8, 10, . . . , (18) and optimized polynomial fits. 

The smooth dependence on the charge densities can be 
parametrized in such a way 

bj ej (pj ,aj) = bj ej (p 3 = 0, a 3 - ) + Pj (af> + 

a^il-p^+a^il-pjf) (4.9) 

that the constraints (|4.2I4.4|) are built in explicitely: 



b x = 1.0, ei(pi = 0,ai) = 0, 
& 2 = l/2, e 2 (p 2 = 0,a 2 ) = -a' 



(4.10) 



In Table UTI we list the coefficients of the fit (I4.9|) . 

In Fig. 0]we compare the first order predictions 12.26fl . 
(|2.33|) for the ground state energies per rung with a Lanc- 
zos diagonalization on a two leg ladder with 8 rungs and 
couplings a = 0.5, a' = 2.7, t' = 0.1(a), 0.2(6), 0.3(c). 
The charge density plateau at p = 1/2 is clearly visible 
in the discontinuous change of the slope in the energy per 
rung as function of p. 
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j 


(0) 




(2) 

°r 


1 


-0.9346 


-1.3769 


0.0167 


2 


0.4154 


-0.2384 


0.0580 



TABLE II: Parameters af A,2) , j 



1,2 in the fit (g^ for the 
energies bjij(pj,ctj) shown in Fig. [3] 



> 

II 





-0.2 
-0.4 
-0.6 
-0.8 
-1 
-1.2 
-1.4 





-•- to- 


(=0.1 




-o- (ft): 


(=0.2 




-A- fcj; 


(=03 










, - ' . -o-^>~~ 










* to 



0.5 
P 



FIG. 4: Zero (dotted line) and first order (solid lines) pre- 
dictions for two leg ladders with a = 0.5, a' = 2.7 and leg 
couplings t' — 0.1(a), 0.2(6), 0.3(c) as well as the correspond- 
ing Lanczos energies for a 2 x 8-site ladder (dashed lines). 



in the plane (2 — a) ft' versus a': 
2- a 



f 



A(a') 

~£ 2 (p2 = 1, a') + £i(pi = 1, a'/2) . (4.11) 



We have determined the chemical potentials pi(pi = 
l,a'/2), p 2 (p2 = from a numerical calculation of 

the ground state energies per site l|4.5(l on systems with 
N r = 8,10,12,14,16,(18) and 

a = 0.0, 0.3,..., 2.7. 

The a'-dependence of A(o/) [right-hand side of (|4.11|) ] is 
shown in Fig. [3J There is a monotonic finite-size depen- 
dence and the thermodynamical limit (solid curve) is es- 
timated with the Bulirsch-Stoer (BST) algorithm. 17 The 
gapped phase with a nonvanishing plateau at p = 1/2 is 
characterized by 



a 



> A(a') 



(4.12) 



i.e. the formation of a plateau at p = 1/2 is favoured for 

• small rung couplings a = J/t < 2 

• large leg couplings a' = J'/t' 



A remarkable agreement between Lanczos diagonaliza- 
tion (dashed curve) and the perturbative result includ- 
ing first order corrections (solid curve) is achieved for 
p < 1/2. A comparison with the zeroth order result 
(dotted curve) demonstrates that first order corrections 
are significant. 

The situation for p > 1/2 is different. Here, we observe 
deviations between the Lanczos diagonalization (dashed 
curve) and the perturbative results (solid curve) which 
increase monotonically with t' and p. Note, that the spin 
coupling J' along the legs (J' = t'a') is already close to 
J' = 1 for t' = 0.3 and a' ~ 2.7. This means in particular 
for p = 1, where the t — J model reduces to a Heisenberg 
model with spin couplings J = 0.5 and J' = 0.81 (for t' = 
0.3 and a' — 2.7), that the product ansatz (|2.14|l with 
rung cluster states (|2.ti|) is an inadequate starting point 
for a perturbative expansion. We expect improvements, 
if we start with a direct product of more complex clusters 
- e.g. plaquettes on the two leg ladder. Moreover, bound 
states (hole-pairs) may change the properties of ground 
state at large charge density (p — > 1). This is indeed the 
case for a > 2ii& 

Let us next turn to the phase diagram of the two leg 
ladder, which is defined by a vanishing plateau width 
W(a,t' ,a') [ H2.37JI ], In general one has to discuss the 
phase boundary in the three dimensional parameter space 
(a, i', a'). In first order perturbation theory in t' [cf. eq. 
(|2.37|) ]. however, it is sufficient to discuss the boundary 




FIG. 5: a'-dependence (first order result) of the plateau width 
of the chemical potential of a two leg ladder at p = 1/2 - finite 
system results and BST-evaluation of the thermodynamical 
limit (TDL). 



B. Three leg ladders 

The ground state energies Ej(pj,aj) of the effective 
Hamiltonians [5th row of Tabled Eas. |TT^|) - (j3~T^|) ] on 
a chain with N r sites have been computed for N r = 
8, 10,..., 16, (18) and 

a = 0.5, a = 2.7 (4.13) 
ai = 1.0125, 0*2 =2.041, a 3 = 6.132. 
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The finite-size dependence has been analyzed with an 
ansatz of the type l|4.5[) for the ground state energies per 
site ay 

The constraints at pj — and at pj — 1 are taken into 
account in the linear interpolations 



-l + 0A Pj 



1,2. 



N 3 (p 3 ) = -1 + 0.232/93- 
The factors bj in Fig. (US) 

h = % 
3 t' 



(4.14) 



(4.15) 



reflect the "renormalization" of the hopping term (6-th 
row in Table QJ. A plot of bjej(pj,atj), j = 1,2, 3, a = 
0.5, a' = 2.7 on the finite systems is shown in Fig. 



a 

CL 



j=2 "' 

OOQ-a--A-0"i?—n O Q-'A-'Q'-" 



v N=8 
N=W 

a N=12 

□ N=14 
o N=16 
a N=IS 



J=3 



Pj 



FIG. 6: Ground state energies per site bji.j(pj,aj) J4.5I4.9|I 
for the effective Hamiltonians Hj (j = 1, 2, 3) of the three leg 
ladder with a = 0.5, a' = 2.7 - results for N r = 8, 10, ... , (18) 
and optimized polynomial fits. 

The data points do not scatter with the system size 
but follow unique curves, which can be considered as a 
reliable estimate of the thermodynamical limit. Their 
dependence on the effective charge densities p e g (fourth 
row in Table^) is parametrized by an ansatz of type l|4.9|l . 
The resulting coefficients with 



h = 1.0 

b 2 = 0.49618 



ei(pi = 0, ai) = 0, 
e 2 {p2 = 0,a 2 ) 



2b 2 



r 2 ^2tb 



b 3 = 0.44028 , e 3 (p 3 = 0,a 3 ) = -a 3 - 



+ 2) 2 
26 2 + 3 



2(6 + 2) 



(4.16) 



are listed in Table ITTll 

Looking at the phase diagram of the three leg ladder we 
have to distinguish four regimes in the three dimensional 
parameter space of a, t', a': 

I plateaux at p = 1/3, 2/3 



j 


«< 0) 




(2) 


1 


-0.7010 


-1.6209 


-0.0260 


2 


0.2054 


-0.4364 


0.1630 


3 


-2.3129 


-0.4460 


0.5739 



TABLE III: Parameters af' h2) , j = 1,2,3 in the fits for the 
energies bjij(pj , ctj) shown in Fig. |S| 



II plateau at p = 1/3, no plateau at p = 2/3 

III plateau at p = 2/3, no plateau at p = 1/3 

IV no plateaux. 

The phase boundaries are defined by the vanishing of the 
plateau width (cf. rows 3 and 8 in Table P): 



P=l/3: 

Wi/3 = £2 - 2ei 



(4.17) 



t' {b 2 fi2{p2 = 1,0-2) + bipi(pi = l,oi)) 



9 = 2/3 : 
W2/3 



e 3 + £1 - 2e 2 



(4.18) 



t' (b 3 p 3 (p 3 = 0, a 3 ) + b 2 p,2(p2 = 0, a 2 )) 



In first order perturbation theory the widths are linear 
in t' and it is therefore convenient to represent the phase 
boundaries in the form 



*i/3 = *i/3 



£2 - 2ei 



2/3 



(bifii(pi = 1, ai) + b 2 p,2{p2 = 1, 02)) 
^ /3 (a,a') 



£3 + £1 - 2£ 2 



(4.19) 



(4.20) 



{b 3 p 3 (p 3 = 0, a 3 ) + 62^2(^2 = 0, a 2 )) 
The four phases defined above are characterized by 
I t f 1/3 (a,a f ) > H > 0, (a, a') > t' > 0, 

II f' 1/3 (a,a') > t' > 0, f > t' 2/3 (a,a'), 

III t' > t' 1/3 (a,a') , t' 2/3 (a,a') > t' > 0, 

IV i' > t' 1/3 (a,a') , t' >t' 2/3 (a,a'). 

We have computed the phase boundaries ^, 3 (a, a') 

and i2/3( a i a ') horn the ground state energies of 

the effective Hamiltonians Hj l|3.12|l on systems with 
10, 12, ... , (18) sites for 

a =1/2 , a' = 0.0, 0.3,..., 2. 7. 

This projection of the phase diagram is shown in Fig. 
El The thermodynamical limits (TDLs; solid curves) are 
estimated with the BST algorithm^ 
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FIG. 7: First order results for the (t', a')-dependence of the 
boundaries of the plateau regimes (I— IV) of a three leg ladder 

- finite system results and BST-evaluation of the TDL. 

V. DISCUSSION AND PERSPECTIVES 

Quasi onedimensional quantum systems - like the 
Heisenberg model or the t — J model - defined on lad- 
ders with I legs show a characteristic sequence of gaps 
in their spectrum. They appear as plateaux in the mag- 
netization curve M(B) and the charge density p{p) (as 
functions of the magnetic field B and chemical potential 
p, respectively) for the spin and charge degrees of free- 
dom. The quantization rule of Oshikawa, Yamanaka and 
Affleck^ defines the values of M and p where these might 
occur. The mechanism which creates the plateaux can be 
studied in a perturbation expansion in the leg couplings 
t'; a' = J' /t' fixed. To zeroth order the ground states 

- at fixed magnetization M and/or charge density p - 
are products of rung cluster states. The latter can be 
classified by a cluster spin £3 and charge Q. First or- 
der perturbation theory leads to an effective interaction 
between the rung cluster states. 

We have studied in this paper the Oth and 1st order 
perturbation expansion in the leg coupling t' (for a 1 fixed) 
on two and three leg ladders. The magnetization - given 
by the total spin - has been assumed to be zero. Our re- 
sults on the two leg ladder can be summarized as follows: 

i) The ground state in the two regimes 

0<p<l/2 (Q = 0,1) 
l/2<p<l (Q = l,2) 

is built up from a direct product of rung cluster states 
with Q = 0, 1 and Q = 1, 2, respectively. 

Even charge states (Q = 0, 2) have total spin and 
are bosonic, odd charge states (Q = 1, (3 for the 3-leg 
ladder)) have total spin 1/2 and are fermionic. For p = 
0, 1/2, 1 all rung cluster states have the same charge Q = 
0, 1, 2, respectively. 

ii) The effective Hamiltonian, which describes the in- 
teraction between the rung cluster states in first order 
perturbation theory, looks like a generalized t — J model 



on a chain, if we treat fermionic states (Q = 1) as elec- 
trons, bosonic ones (Q = 0, 2) as holes. A comparison 
of the perturbative results with exact diagonalizations 
yields good agreement for the ground state energies in 
the first regime (p < 1/2, t' < 0.3), but increasing devi- 
ations with p and t' in the second one. 

Hi) Increasing the size of the cluster we started with, 
will improve these results and might reveal further gaps 
in the spectrum. In a next step we will consider a direct 
product of plaquette cluster states on the two leg lad- 
der (Fig. P). Here, cluster ground states with charges 
Q(x) = 0,1,2,3,4 occur, which might induce additional 
plateaux at p = 1/4 and p = 3/4. Indeed, evidence for 
the existence of a long range charge density wave state 
in the t — J two leg ladder at p = 3/4 has been found in 
a DMRG calculation^ 

iv) We have studied the stability of the charge density 
plateau and the phase diagram in the regime where first 
order perturbation theory is applicable. We found that 
the formation of a charge density plateau is favoured for 

• small values of the couplings a — J/t on the rungs 

• large values of the couplings a 1 = J' jt' on the legs. 

In the three leg ladder case, the ground states in Oth 
order are built up again from rung cluster states with two 
charges: 

0<p<l/3 (Q = 0,1) 
l/3<p<2/3 (Q = l,2) 
2/3<p<l (Q = 2,3) 

For p = 0, 1/3,2/3, 1 all the rung clusters have the same 
charges: Q = 0, Q = 1, Q = 2, Q = 3, respectively. 

The method developped here for ladder systems should 
be applicable whenever finite clusters containing the 
"large" couplings can be defined in a natural way. If 
we for instance consider the t — J model on a Shastry- 
Sutherland lattice^ the "large" couplings on the diag- 
onals define two site clusters. The same arguments we 
developped here for the two leg ladder (with two site 
clusters) can be repeated and we expect again a charge 
density plateau at p = 1/2. Of course the effective inter- 
action between the cluster states depend on the geometry 
of the lattice and are therefore different for the Shastry- 
Sutherland lattice and the two leg ladder. 

The method developped here for charge density 
plateaux is also applicable to magnetization plateaux 
provided that finite clusters containing the "large" cou- 
plings can be defined in a natural way. Again, the Heisen- 
berg model on a Shastry-Sutherland lattice^ is a good 
example. 

However, for a realistic description of the experi- 
mentally found magnetization plateaux M — S/N = 
1/6, 1/8, 1/16 (S* total spin, N total number of sites) one 
should start from clusters which at least - requiring inte- 
ger total spin for the clusters - contain 6, 8, 16 sites! Of 
course the computation of the cluster ground states and 
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in particular of their interaction becomes more involved 
for decreasing magnetization values. 
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